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Abstract: We theoretically study macroscopic entanglement between 
a magnetically trapped Bose-Einstein condensate and a superconducting 
loop. We treat the superconducting loop in a quantum superposition of two 
different flux states coupling with the magnetic trap to generate macro- 
scopic entanglement. The scheme also provides a platform to investigate 
interferometry with an entangled Bose Einstein condensate and to explore 
physics at the quantum-classical interface. 
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1. Introduction 

Entanglement is considered to be one of the most fundamental features of quantum mechan- 
ics. In addition, it is of great importance in the context of quantum information and quantum 
computation. In recent years, there have been considerable efforts to generate and preserve 
entanglement for quantum information processing. In particular, entanglement among macro- 
scopic observables is of prime interest in order to explore how a system behaves at the interface 
of classical and quantum mechanics. Microscopically entangled states are also promising can- 
didates for the practical realization of a quantum computer. 

In recent years there has been ground breaking progress in the field of manipulation of BECs. 
Nowadays, it is relatively easy to produce a BEC in a micro-magnetic trap on an atom chip (T|. 
A BEC in such traps can be coherently manipulated with RE fields J2] [3] and microwaves 
HE). Since neutral atoms can be positioned a few microns from the chip surface and moved 
with nanometre resolution, atom chips provide a convenient platform to study the interaction 
between a BEC and a nearby surface |6) including the study of fundamental quantum effects 
such as the Casimir-Polder interactions |7 j. On the other hand the field of superconducting cir- 
cuits is progressing rapidly in terms of technological implementation and realization of quan- 
tum coherent control of superconducting qubits EH). A macroscopic superposition of different 
magnetic flux states has been demonstrated iflOl [TT1 and quantum coherent dynamics of flux 
qubits have been realised (12) . In addition, there are proposals to study the interaction between 
a superconducting circuit and a micromechanical device, in particular to entangle a microme- 
chanical resonator with a Cooper pair box |[T3l . Since a nanomechanical cantilever can be easily 
integrated on an atom chip a magnetically trapped Bose-Einstein condensate can be coherently 
coupled to the oscillation modes of a magnetized cantilever fl4l . Two spatially separated Ry- 
dberg atoms can also be capacitively coupled to each other through a thin superconducting 
wire fT5) ■ Atom chips made of superconducting substrate and wires carrying persistent current 
to manipulate trapped atoms have been reported lfl6l[T7ll . The magnetic field exclusion due to 
Meissner effect in a superconducting wire can significantly modify the magnetic trap properties 
J25). Recently, there is a great interest to study the coupling between ultracold gases and solid 
state quantum systems to explore their quantum properties. We theoretically study a macro- 
scopic entanglement between a superconducting loop and a BEC. By exploiting the quantum 
mechanical properties of superconducting circuits and their interactions with nearby ultracold 
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Fig. 1. Schematic showing a superconducting loop situated symmetrically above a Z- 
shaped wire magnetic trap. The persistent current in the loop can exist in a quantum super- 
position in the clockwise and anti-clockwise directions. The magnetic field produced by the 
persistent current can perturb the BEC trapping potential where the degree of perturbation 
can be increased by moving the trap center closer to the loop. The superconducting loop 
can be biased appropriately by applying an external magnetic field along the z-direction. 

atoms or a BEC a macroscopic entanglement can be produced which provides a pathway to 
explore physics at the quantum-classical interface. 

2. Superconducting loop on an atom chip 

According to the present scheme, a superconducting loop existing in a quantum superposition 
of macroscopically distinct flux states can perturb the BEC trapping potential and under certain 
conditions a macroscopic entanglement can be generated. The macroscopic variables involved 
in the entangled state are the flux linked to the loop, the chemical potential of the BEC and the 
position and shape of the BEC wavefunction. The physical arrangement of the superconducting 
loop on the atom chip is shown in FigQ] Ultracold neutral atoms or a BEC can be trapped in 
a magnetic trap above the surface by applying an external bias field in the plane of a Z-shaped 
current-carrying wire. A superconducting loop is positioned symmetrically above the Z-wire 
and below the trap in such a way that there is zero net flux linked to it. The trap position and trap 
frequencies depend on the external bias fields and the current in the Z-wire. The magnetic field 
produced by a nearby superconducting loop carrying persistent current can perturb the magnetic 
trap potential through magnetostatic interactions. The sign and amplitude of the perturbation 
depends on the direction of the persistent current flowing in the loop and its position with 
respect to the trap centre. 

A superconducting loop when placed in an external flux permits only discrete values of the 
net flux threaded through it, which is an integral multiple of the flux quantum. In other words, 
the super-current in the loop responds automatically to any change in the externally applied 
flux in order to keep the closed loop phase acquired by the wave function an integral multiple 
of 2k . In the case of a superconducting ring interrupted by a Josephson tunnel junction the 
total energy (the sum of the magnetic energy of the loop and Josephson coupling energy of 
the junction) corresponds to a double-well potential of a flux variable threaded through the 
loop ITOl . where the left (right) well corresponds to persistent current flowing in a clockwise 



(anticlockwise) direction in the loop. The double well becomes symmetric when the applied 
flux is equal to half of the flux quantum. The barrier height between the wells depends on the 
critical current of the Josephson tunnel junction. The inter-well tunnelling can be controlled 
externally by replacing the single Josephson junction with split junctions (DC-SQUID) flTOl . 

At low temperature the superconducting loop can be prepared in a quantum superposition 
of two persistent current states flowing clockwise and anticlockwise by biasing it at half of 
the flux quantum. A BEC can be prepared in a Z-shaped wire magnetic trap far above the 
superconducting loop where the perturbation in the trapping potential caused by the loop is 
negligibly small. If the magnetic trap containing a BEC is moved adiabatically towards the 
superconducting loop carrying a persistent current the trapped BEC will adiabatically follow the 
perturbation in the potential induced by the magnetic field from the loop. However, if the loop is 
initially prepared in a quantum superposition of persistent current flowing clockwise and anti- 
clockwise it will perturb the magnetic trap potential in two different perturbations which can be 
distinguished by increasing the interaction with the loop leading to a macroscopic entanglement 
between the persistent current state of the loop and the state of the BEC in different perturbed 
potentials in the magnetic trap. One perturbed configuration of a BEC can differ from the other 
one in terms of its spatial distribution and the chemical potential. The signature of entanglement 
between the superconducting loop and the BEC can be observed through the atomic distribution 
in time of flight measurements. 

3. Entanglement between the loop and a BEC 

In order to study the macroscopic entanglement between the loop and the BEC we start by 
deriving the Hamiltonian of the superconducting loop coupled to the magnetic trap containing 
a BEC. In the case when the loop is biased at half of the flux quantum the double well governing 
the dynamics of the loop is symmetric. Therefore, at low temperature, the Hamiltonian of the 
superconducting loop can be treated as a two-level system. 

ff s = £o|0)(0|+£ |l)(l|+/|l)(0|+/|0)(l| (1) 

where |0) and 1 1) represent the ground state of the left and the right well, respectively, and J is 
the tunnelling amplitude between them. 

The Hamiltonian of an atom of mass m in the ground state of the trap coupled to the super- 
conducting loop in the case 7 = can be written as 

H T = /V(r) 

where AVb(r,f) and AVi(r,t) are the perturbations in the trap potential V(r) due to coupling 
to the loop in states |0), 1 1), respectively, and |0) (0| + 1 1) (1 1 = /. For a magnetic trap potential 
V(r) = mfgfHgB(r) the perturbation in V(r) originates from the magnetic field produced by 
the persistent current in the loop which affects the harmonic trapping field profile B(r) of the 
Z-wire magnetic trap near its trap centre. In the case when loop is prepared in state |0) and 
the perturbation is increased adiabatically the trapped BEC follows the ground state ^Q(r,t) 
of the trapping potential V(r) + AVo (r,f) and likewise if the loop is initially prepared in state 
|1) the BEC follows the ground state 0i(r,f) of the potential V(r) +AVi(r,t). However, under 
the adiabatic condition and for the loop in quantum superposition (of |0) and |1)) the trapped 
BEC follows two different configurations <poi r ,t) and 0i (r,t ) because the state of perturbation 
is coupled to the state of the loop. Therefore, the field operator ^(r,?) can be expanded as a 
linear combination of <po(r,t)\0) (0| and 0i(r,f)|l)(l|; 



— V z + V(r) )/ + AVb(r,f)|0)(0|+AVi(r,r)|l)(l| 



V>(r)dr (2) 



❖(r,O=flo^(r,O|O)(O|+fli0i(r ) r)|l)(l| (3) 

where &q and a\ are the corresponding bosonic annihilation operators. Therefore, from 
Eqs.[T][2]and[3]the Hamiltonian of the system can be written as 



where 



and 



//(r)=£ |0)(0|+£o|l)(l|+Aio(f)4«o|0)(0|+ J Lii(0«1«i|l)(l| (4) 
MO - / 0oV) (^ v2 +V(r) + AV (r,t)\ 0o(r,f)dr (5) 
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(6) 



are the energy eigenvalues in the case of two perturbed situations of the trap. 

Let us see how entanglement can be generated. In the first step a BEC of N atoms is prepared 
in a Z-wire magnetic trap far away from the superconducting loop so that perturbation caused 
by the loop is negligibly small. In the second step the superconducting loop is prepared in a 
symmetric superposition \S) = (|0) + |l))/v2 and then the tunnelling amplitude J is reduced 
to zero. In the third step the BEC is slowly brought closer to the superconducting loop (by 
increasing the jc-bias magnetic field or by decreasing the current in the Z-wire) so that the per- 
turbation to the magnetic trap is increased adiabatically and the BEC will follow two different 
configurations of the trap. As the two different perturbations grow in quantum superposition 
a macroscopic entanglement between the state of the BEC in the trap and the state of the su- 
perconducting loops established. At one point where the coupling is sufficiently strong two 
perturbed configurations of the trap will be different and distinguishable. 

Therefore at t = the initial state of the system is \m,t = 0) = (|0) + 1 \))\N, 0(r,f =0))/V2, 
where \N, (j)(r,t = 0)) is the state corresponding to N atoms in the ground state <j>(r,t) of the trap 
potential V(r) in the case of no coupling with the flux loop. When the coupling is increased adia- 
batically the state of the system evolves under the unitary evolution U (f) = exp(— i J H(t')6t') 
to 



„iy a (t)-i(E a t+N S< ^(t')di')/h 



V2 

where 



|0)|Ar,^(r, f ))+ e i *«|l)|AT^ 1 (r,/)) 



(7) 



rf ^^W ttW-ttft) (8) 







and \N,(po(r,t)) and \N,<j>i(r,t)) are states corresponding to N atoms in the ground states 
0o(r,f) and 0i(r,f), respectively. In addition 70(f) and 71(f) represent the geometrical phase. 
Equation. |7]is derived using the following transformations implied by the adiabatic condition 

e -|/oW)('')«5 d bl )(°l d ''|o)|Ar,(|)(r,/ = 0)) =«-^^W)(Od«' e «»(0|o)|^,^ ) (r,0> (9) 

e -^Pi(O^i|i)(i|^|i)|^ > ^( r> ( = 0))=e-¥^/ t i(»0d«' e «n(0|i)|^ > ^ 1 ( r ^)) (10) 

It is evident that l^f) (Eq .|7]i represents a macroscopic entangled state between the states 
of the loop and the BEC in two different configurations. As <t>(f) evolves l^f) spans different 
macroscopically entangled states. 
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Fig. 2. Magnetic field intensity profile along the axial direction of a Z-wire magnetic trap 
coupled to a superconducting loop when the persistent current flows (a) clockwise and (b) 
anti-clockwise. The dotted curve (red) represents the harmonic field profile without any 
coupling. The dashed curve (green) and solid line (blue) show the effect of perturbation 
when a persistent current in the loop generates a flux through the loop of magnitude one- 
quarter and one-half of a flux quantum, respectively. . 



In order to calculate the perturbation in the magnetic field intensity profile of the Z-trap a 
circular superconducting loop of diameter 10 jJ.m is considered. The loop is positioned at a dis- 
tance 10 jiim below the trap minimum with its axis parallel to the z-axis. In order to bias the loop 
to produce a symmetric energy double well an applied field should generate a flux of magnitude 
one-half of a flux quantum through the loop. For the loop parameters given above a uniform 
magnetic field of about 0.1 G is required in the z-direction to bias the loop. The trap parameters 
are calculated for a Z-wire magnetic trap carrying 5 A DC current and 20 G uniform trapping 
bias magnetic field in the x-direction. The trap bottom is tuned to 1 G by applying an additional 
uniform bias field along the y-direction. The perturbation in the magnetic field intensity profile 
of a Z-wire trap in the axial direction due to the magnetic field from the superconducting loop 
is shown in Fig. |2] for clockwise (a) and anticlockwise (b) flow of persistent current (same in 
magnitude). The dotted curve (red) represents the unperturbed harmonic profile when the per- 
sistent current is zero. The dashed line (green) and solid line (blue) represent the magnetic field 
intensity corresponding to the persistent current generating a flux through the loop of magni- 
tude one-quarter and one-half of a flux quantum, respectively. A three-dimensional profile of 
the magnetic field intensity (around the trap centre) in a horizontal plane parallel to the x-y 
plane and passing through the trap minimum is shown in Fig. [3] The axial and radial axes refer 
to the weakly and strongly confining directions, respectively, in the horizontal plane. The cor- 
responding density plot is shown in Fig. [4] The amplitude of the perturbation which is defined 
as the difference of the magnetic field intensity at the local maximum and the local minimum is 
about 5.5 mG for a persistent current generating a flux of magnitude one-half of a flux quantum 
through the loop as indicated in Fig. 12(b). The chemical potential for N atoms of 87 Rb (for 
\F = 2,m,F = 2) with radial and axial trap frequencies 540 Hz and 10 Hz, respectively) in the 
absence of any perturbation is 0.01959./V 2/ ' 5 mG. The amplitude of the perturbation should be 
greater than the chemical potential of the BEC in order to distinguish between the perturbed 
states. The distance between the minima of the two perturbed magnetic field profiles is of the 
order of the diameter of the superconducting loop. The strength of the perturbation depends on 
the distance of the magnetic trap centre from the superconducting loop and increases as the trap 
centre is displaced towards the loop. Figure [5] shows the perturbation amplitude as a function 
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Fig. 3. Magnetic field intensity profile around the trap minimum when a clockwise flow of 
persistent current generates a flux through the loop of magnitude one-half of a flux quan- 
tum. Radial and axial axes lie in the horizontal plane passing through the trap minimum. 
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Fig. 4. Density plot showing the magnetic field intensity in the horizontal plane. The colour 
scheme has one to one correspondence with that indicated in Fig. [3] 
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Fig. 5. Perturbation amplitude as a function of distance (d) between the superconducting 
loop and the centre of the magnetic trap when the persistent current generates a flux through 
the loop of magnitude one-half of a flux quantum. 



of distance (d) between the trap centre and the loop measured along the axis of the loop. It 
is evident from Fig. [5]that the perturbation amplitude at d > 20 jim is negligibility small and 
increases rapidly as d is reduced. 

4. Entangled atom interferometry and detection 

The interference property of a Bose Einstein condensate is an important tool to study the co- 
herence among the condensates. In an analogue of a double slit experiment a single BEC can 
be coherently split by gradually changing the harmonic confinement into a double well and an 
interference pattern can be observed by overlapping two split BEC wavefunctions after they 
are released ll23l 0. In such experiments BECs are not path entangled. Now it is important 
question: if we release the BEC from macroscopically entangled state (Eq. |7]i can we get an 
interference and how can we detect macroscopic entanglement between the superconducting 
loop and the BEC? 

In order to make the idea clear let us consider a single atom in the BEC (N = 1). The expec- 
tation value of the atomic density after the release is evaluated as 

«(r,f) = (^,f|^t( r )v/( r )|xp, f ) = I(|0 o (r,O| 2 + |^i(r,r)| 2 ) (11) 

This is an incoherent mixture with no interference; the interference term is suppressed be- 
cause the which-path information about the BEC state is available through the entanglement 
with the superconducting loop. It is possible to know which trap configuration the atom oc- 
cupies just by measuring the state of the loop. Therefore, in order to get interference the 
which-path information should be erased. This is possible if |0) and |1) are indistinguish- 
able. On the other hand, in a new basis where |+) = -^(|0) + |1)), |-) = -^(|0) - 

\S) = ±(\N=lMr,t)} + \N = l,<h(r,t))) and |A) = ±(\N = l,fc(r,t)) -\N= 
the state Eq.|7]for <t>(f ) = can be written as (overall phase ignored) 



\^,t) = ±[\+)\S) + \-)\A)} (12) 

It is evident from the state Eq. [12] that when the superconducting loop is found in state |+) 
in repeated measurements the corresponding atomic distribution shows an interference term 
cos(27Tr/A), where 2n/A = tmd/(h(t 2 + (mOg /h) 2 )) (initial gaussian width Oq located at 
id/2) ||23l . If the loop is found in | — ) the atomic distribution exhibits an anti-interference term. 
Therefore, if the state of the superconducting loop is not measured then there will be no way to 
differentiate between interference and anti-interference and no interference pattern will emerge. 
In the case of more than one atom in the state Eq. [JJ the situation is different. If we treat the 
BEC state as the control qubit acting on the superconducting state and perform a CNOT trans- 
formation (UcNor\N,Mr,t))\0) = \N,h(r,t))\0) and Uatar\N,h(r,t))\l) = \N,h(r,t))\0)) 
after releasing the BEC then the state Eq.[7jtransforms to 



\V,t) = -J=|0)(|^V^o(r,f))+ e 1 '*(' ) |^V^ 1 (r,0» (13) 

There is no entanglement between the BEC and the superconducting loop now and in this 
case the centre of mass density of the BEC when released from the state Eq. Qj] should exhibit 
a pattern with A = A/N lfl"8l [191 l20l |2"T1 l22l in repeated time of flight measurements; this is 
a signature of the macroscopic entanglement Eq. [7j However, only those measurements which 
count exactly the same N should be considered; therefore, a precise control for state preparation 
(Eq.[7j is required, where <&(/) should not vary by more than over n/2N from measurement to 
measurement l22l . However, it might be difficult at present to implement such a macroscopic 
CNOT operation in order to get the state of interest Eq. Q~3] There is another way to observe 
a signature of the macroscopic entanglement between the loop and the BEC. If we measure 
the distribution of the BEC when it is released from the trap and correlate it with the state of 
the loop then it is evident from Eq. |7] that whenever the superconducting loop is measured in 
|0) the measured particle density distribution should be N\<f>o(r,t)\ 2 and when it is measured 
in 1 1 ) the particle density distribution is given as N\fa{r,t)\ 2 . This can be evidence that the 
superconducting loop is interacting with the BEC in the trap. However, if the superconducting 
loop is measured in state |+) then in this basis for the state Eq.|7Jthe particle density distribution 
should shift either from N\(j>o(r,t)\ 2 or from N\(j>i (r,t )\ 2 J22). This shift in density distribution 
can be attributed to the existence of macroscopic entanglement of a pure condensate with a 
superconducting loop. 

5. Decoherence 

In context of a practical realization it is important to mention that the generation of macroscopic 
entanglement demands the validity of adiabatic turn-on of the perturbation potential. Therefore, 
the decoherence time for the loop must be larger than the time required to complete the adiabatic 
process so that the required entanglement can be observed. There are various factors 11241 which 
can limit the decoherence time of a flux loop. Noise in the magnetic field, interaction of the 
loop with the state measuring apparatus and impurities in the substrate on which the loop is 
fabricated can decohere a flux superposition. At present the decoherence and relaxation times 
for a flux qubit of the order of 1-10 /is have been reported ll26l . In order to satisfy the adiabatic 
condition while increasing the coupling between the superconducting loop and the magnetic 
trap the turn on process should be slow enough to satisfy ^ <C CO 2 where CO is the angular trap 
frequency. If change in the trap frequency during the perturbation turn-on process is of the order 
of the trap frequency then the time for the adiabatic approximation to hold can be estimated 



roughly from the trap frequencies. Since the effect of perturbation is much stronger along the 
axial direction (weaker confinement). Therefore, the inverse of axial trap frequency can give 
a rough estimation of the adiabatic time limit. For a typical axial trap frequency of 10 Hz the 
adiabatic time limit estimation is about 100 ms which requires a considerable improvement in 
the decoherence time of the flux loop. However, an accurate calculation of the time limit for the 
validity of the adiabatic approximation can be calculated by numerically solving the dynamics 
of BEC in the time dependent perturbation by minimizing the transition probability to other 
excited states from the instantaneous ground state of the potential. 

In order to evaluate the time limit for the adiabatic approximation to hold it is important 
to know the functional form of the magnetic trap potential moving towards (z-direction) the 
superconducting loop which can be expressed as 

V(x,y,z,t) = ^m(0 2 (t)x 2 + ^m(0 2 (t)y 2 + ^m(0 2 (t)(z + zo(t)) 2 + V + 

^-£|0><0|-^-$|l)<l| (14) 

where co x , (O z correspond to the radial trap frequencies and co y is the axial trap frequency and it 
is assumed that the state dependent perturbation is well defined. Since the trap considered here 
is predominantly anisotropic (co Xl co z >> (fly) the effect of perturbation induced by the loop is 

more significant along the axial direction. The term 2a ^} y e a 7 represents the mean perturbation 
along the axial direction of the magnetic trap when the loop is in 1 1) state. When the loop exists 
in state |0), the sign of the perturbation is changed. The displacement of the trap towards the 
loop and how fast the perturbation is induced depends on zo(t). The trap frequency in the case of 
a Z-wire trap is inversely proportional to the distance of the trap centre from the Z-wire and the 
offset field at the trap centre. A variation in the trap frequency when the trap is moved towards 
the loop can be avoided by synchronously varying the offset field. According to the numerical 
calculations shown in Fig.[5]the perturbation amplitude is negligibly small for d > 20 /im and 
grows rapidly as d is reduced. Therefore a BEC can be prepared at about d > 20 jim and by 
decreasing d to about 10 jim a perturbation amplitude of about 5.5 mG can be obtained. The 
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function Bq + koy 2 + a o is fitted to the axial magnetic field profile shown in Fig. 12(a) 

°o 

for clockwise flow of supercurrent generating a flux equal to half of a flux quantum (solid 
blue line) where values of the fit parameters are Bo = 999.85 mG, ko = 0.00031 mG//im 2 , 
(To = 10.13 jJ.m and a = —32.0 mGfim. Therefore, the potential profile can be obtained by 
using V(r) = mfgFjiBB(r). The sign of a changes for the anti-clockwise flow of persistent 
current. One can use the potential described in Eq.[14]in order to study the dynamics of BEC 
coupling to the superconducting loop and to numerically evaluate the time limit for the adiabatic 
approximation. 

6. Remarks 

In the context of experimental realization it is important to consider that the magnetic field 
from the Z-wire and the bias field should be less than the critical field of the superconducting 
loop. The experiment can be constructed by utilizing flip-chip technology, where the Z-wire 
and the superconducting loop can be constructed on two different substrates which can be 
bonded together with high precision with an appropriate gap between them. The Z-wire can 
also be constructed from a superconducting material in order to reduce the technical noise in 
the current. The atom chip can be shielded from the background radiation by a gold coated 
copper shield. The whole assembly can be mounted on a cold finger with the chip pointing 



up side down. The ultracold atoms can be prepared in a different chamber and magnetically 
transported to radiation shielded chamber containing the chip where they can be trapped and 
evaporatively cooled ifTTI down to a BEC. 

There are various factors that can destroy a macroscopic entangled state described by Eq. [7] 
Loss of a single atom from this state can easily destroy the entanglement by sharing its infor- 
mation with the environment. However, a superconducting flux superposition itself is prone to 
environment-induced decoherence (8). In order to satisfy the adiabatic condition while increas- 
ing the coupling between the superconducting loop and the magnetic trap the coupling turn-on 
process should be slow enough but should be completed before the decoherence time limit of 
the superconducting loop otherwise excitations can occur in the BEC. 

In conclusion we have shown how the coherent dynamics of a superconducting loop can be 
used to generate a macroscopic entanglement of a BEC on an atom chip. Such a macroscopic 
entanglement could be useful to explore fundamental quantum mechanics by studying how 
quantum mechanical effects behave at the macroscopic level and how they decohere with the 
system parameters. It may also be possible to explore decoherence between the superconducting 
circuit and the BEC by varying the trap parameters and the size of the system such as the number 
of atoms in the BEC. It has been shown how to realize interference of an entangled BEC and 
in addition how it can be used to detect a signature of macroscopic entanglement between the 
superconducting loop and the BEC. 
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